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We present results on spectra of pions, kaons and (anti)protons from a study of heavy 
ion collisions using the perturbative QCD + saturation model to calculate the production 
of initial (transverse) energy and baryon number followed by a hydrodynamic description 
of the expansion of produced matter. In particular, we study how the hadron spectra 
and multiplicities depend on the decoupling temperature T& ec when the low temperature 
phase contains all hadrons and and hadron resonances with mass below 2 GeV. We show 
that the spectra and multiplicities of pions, kaons and (anti)protons measured at RHIC 
in central Au+Au collisions with \fs = 130 GeV can be obtained with a single decoupling 
temperature 150. . .160 MeV, common for both the chemical and the kinetic freeze-out. 



1. INTRODUCTION 

The average transverse momentum of primary particles in heavy ion collisions is ex- 
pected to be larger than the experimentally observed one. Two ways out of this problem 
have been discussed: In one case the number of produced particles is small and the fi- 
nal multiplicity is achieved after fragmentation which reduces also the average transverse 
momenta of final, observed particles [ |Ij]. On the other hand, results from the calculation 
can be used as initial conditions in a hydrodynamical approach [ ||. In this case the extra 
transverse energy is transferred into the longitudinal motion through the work by pressure 
in the expansion. The transfer from transverse to longitudinal degrees of freedom takes 
place because the matter is produced in a strong longitudinal expansion. 

When time passes, also a transverse flow builds up and leads, e.g., to a clear mass 
dependence of final hadron spectra. The shape of spectra depend both on the transverse 
flow and the decoupling temperature Td ec but the ratios of particle multiplicities on the 
Tdec and on the baryon chemical potential /i_B,dec- Strangeness chemical potential is fixed 
by demanding that the net strangeness density is zero. 

Modelling the heavy ion collisions at collider energies in terms of perturbative QCD 
calculation of minijet production at saturation [ [| and of hydrodynamic expansion of 
produced matter, has lead to a successful description of multiplicities and transverse 
energy at mid-rapidities [ [3[ . This talk summarizes the predictions for particle spectra in 
central Au+Au collisions at y/s = 130 GeV from the pQCD+saturation+hydrodynamics 
model [ |] . 
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2. THE MODEL FOR CALCULATING THE TRANSVERSE SPECTRA 

In a high energy AA collision parton production can be calculated from the nuclear 
parton distribution functions and the perturbative parton-parton cross sections, provided 
that the momentum transfer in the partonic collision px S> Aqcd- Assuming independent 
partonic collisions, the dominant part of the initial production is obtained by extending 
the (LO) calculation of partonic jets down to po — 1 . . . 2 GeV3> Aq C d [ 0- 

The basic quantities for characterizing the parton production in nuclear collisions are the 
number and the total transverse energy of minijets per unit rapidity. These are obtained as 
a function of p$ from the transverse energy distribution of minijets in the central rapidity 
unit as the normalization cjjet(po) and the first moment a (Ex) [ §. Combining these with 
the nuclear overlap function Taa (or the product of thickness functions of colliding nuclei) 
gives the average total number and transverse energy of partons with p^ > Po- The net 
baryon number can be obtained as a difference of the quark and anti-quark numbers [ [J. 
Finally, to close the model, a cut-off momentum must be fixed and here we assume that it 
originates from the saturation of partons [ [7] . We express this in terms of final, produced 
partons as a simple geometrical condition: each parton is given a transverse size tt/p t , 
and when the partons overlap, extra independent primary collisions with smaller px are 
suppressed. The lower cut-off po = p sa ,t is thus fixed from the saturation condition [ || 
N(p ,Ay)-7c/p 2 = nR 2 A . 

From the minijets at saturation we can calculate the energy density and the net baryon 
number density after introducing a correlation between the longitudinal momentum and 
the longitudinal space-time formation point of the minijet. We assume that rapidity of the 
minijet equals the space-time rapidity of the formation point, y = r\ = 0.5 ln[(t+z) / (t—z)). 
The formation time, defined as l/p sa .t is taken as the initial time for the hydrodynamic 
evolution. The details of the determination of the initial state can be found in [ |], [3| . 

It can be argued that r = l/p sa t (~0.2 fm/c for Au+Au at RHIC) is too short to achieve 
thermalization. However, from the point of view of the average energy per particle, the 
system does look thermal [ 0. Also, one should expect collisions to take place and the 
build-up of collective behaviour to start already during the thermalization process. From 
this point of view, the use of hydrodynamics at early times should be considered to 
approximate the effects of momentum transfer in the collisions during equilibration. 

To solve the hydrodynamic equations, an equation of state (EoS) is needed. We assume 
an ideal QGP high temperature phase with a first order phase transition to a hadron 
resonance gas at T c = 165 MeV. All hadrons and hadron resonances with M < 2 GeV are 
included and a repulsive mean field is assumed in order to ensure a consistent temperature 
behaviour of the pressure [ ||. The final spectra are obtained using the Cooper-Frye 
freeze-out procedure and the two- and three-body decays of all resonances are included. 

3. RESULTS 

In Figure 1 the Td ec dependence of spectra of negative pions, kaons and antiprotons is 
shown together with the data measured by the PHENIX Collaboration [ [|. The slopes of 
the spectra are seen to change considerably with Td cc - For px > 1.5 GeV, the decrease of 
T dec from 150 MeV to 120 MeV changes the inverse slope T efr from 265 MeV to 305 MeV 
for pions, from 275 MeV to 330 MeV for kaons and from 295 MeV to 375 MeV for protons. 
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For the (pt) the changes are from 0.46, 0.65 and 0.86 GeV to 0.53, 0.78 and 1.08 GeV for 
pions, kaons and protons, respectively. The change of the decoupling temperature affects 
also the normalization of spectra, that is the multiplicity. In fact the multiplicity of pions 
increases slightly whereas that of kaons decreases 21 % and of protons 52 %. As has 
been observed earlier [ [n|], to reproduce in thermal models the multiplicities of heavier 
particles and pions simultaneously, a freeze-out temperature for chemical reactions must 
be of the order of 150 to 160 MeV. In our calculation not only the multiplicities but also 
the observed shapes of the spectra are obtained if the kinetic freeze-out and the chemical 
freeze-out are assumed to take place at the (approximately) same temperature of order 
150. . .160 MeV. 




Figure 1. Spectra of negative particles cal- 
culated [g for T dcc = 160, 150 and 120 MeV. 
The data points are from a measurement by 
the PHENIX collaboration f 01. 




Pt [GeV] 

Figure 2. As in Fig. 1 but for positive par- 
ticles and with T^ ec = 150 MeV only. The 
feed-down from hyperon states is included in 
the data but not in the calculation. 



In Figure 2 the spectra of positive pions, kaons and protons are depicted using T dec = 
150 MeV in the calculation. Within the accuracy of the data the normalizations of pions 
and kaons agree but the calculated proton spectrum is slightly below the measured one. 
This is to be expected since the feed-down from hyperons is included among the measured 
but not the calculated protons. PHENIX has also reported p and p yields corrected for 



hyperon feed-down, dN/dy = 19.3 ± 0.6 and 13.7 ± 0.7 [ |TT| . Our results for Td ec = 150 



MeV are 20.2 and 13.1 in excellent agreement with the measurement. Also the ratio 
A/A = 0.70 is consistent with the PHENIX result 0.75 ± 0.09 [ [IJ. 

Note that the deviation of p/p and A/A from unity is due to the net baryon number 
content of the initial matter at saturation. In the pQCD + saturation calculation the 
net baryon number cannot be changed independently of the total multiplicity. The rela- 
tive production of quarks and antiquarks (and gluons) is completely fixed by the parton 
distributions and the perturbative parton level cross sections. In addition to the net 
baryon number, the number of baryons relative to that of antibaryons depends strongly 



4 



on decoupling temperature and again in our calculation agreement is obtained only for 
T dpr ~ 150. . .160 MeV. 



4. DISCUSSION 

Using the pQCD + saturation model [ |J to calculate the initial particle production in 
central Au+Au collisions at \/s = 130 GeV and hydrodynamics to describe the expansion 
of the matter we obtain a good agreement between the calculated and measured quantities. 
We find it remarkable that in order to reproduce either the multiplicities, the slopes 
or the antibaryon-to-baryon ratios we come up in each case with the same decoupling 
temperature of order 150. . .160 MeV. Similar conclusion was also obtained in [ [12| . 



It has been argued that To = l/p sa ,t ~ 0.2 fm/c is too short a time for achieving full 
thermalization. This might well be so but it is not plausible to assume that there are 
no collisions with momentum transfers among the constituents at times before a possible 
later thermalization time. Collective motion, the flow, is a result of the collisions among 
the particles of the matter which lead to a net momentum transfer from denser to less 
dense parts in the matter. From this point of view it is difficult to justify the starting 
of a hydrodynamical calculation at a later time assuming no initial transversal velocity. 
In any case, we have checked the effect of the initial time on the shapes of the spectra. 
Taking r = 0.6 fm/c, needed to explain the formation of elliptic flow [|R|, while keeping 
the total entropy fixed, changes the slopes only slightly. Even r = 1.0 fm/c is not enough 
to compensate for the change from the decrease of Td ec from 150 to 120 MeV. 

Finally, we want to emphasize that our model scheme is very tight. Uncertainties in 
the pQCD + saturation calculation are essentially fixed from a single global quantity, like 
the total multiplicity. In the hydrodynamical part a large number of hadrons and hadron 
resonances must be included but after fixing the EoS, the only real freedom left is the 
decoupling temperature Td cc - 
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